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GENERATING INVISCID AND VISCOUS 
FLUID FLOW SIMULATIONS OVER A 
SURFACE USING A QUASI-SIMULTANEOUS 
TECHNIQUE 

5 

STATEMENT REGARDING FEDERALLY 
SPONSORED RESEARCH OR DEVELOPMENT 

This invention was made with Government support under 
contract NNL08AA08C awarded by NASA. The Govern- to 
ment has certain rights in the invention. 

BACKGROUND 

1. Field 15 

This application relates generally to simulating a fluid flow 

over a surface and, more specifically, to generating both invis- 
cid and viscous fluid-flow simulations using a quasi-simulta- 
neous technique for a fluid flow over a computer-generated 
aircraft surface. 20 

2. Description of the Related Art 

Aerodynamic analysis of an aircraft moving through a fluid 
typically requires an accurate prediction of the properties of 
the fluid surrounding the aircraft. Accurate aerodynamic 
analysis is particularly important when designing aircraft 25 
surfaces, such as the surface of a wing or control surface. 
Typically, the outer surface of a portion of the aircraft, such as 
the surface of a wing, is modeled, either physically or by 
computer model, so that a simulation of the fluid flow can be 
performed and properties of the simulated fluid flow can be 30 
measured. Fluid-flow properties are used to predict the char- 
acteristics of the wing including lift, drag, boundary-layer 
velocity profiles, and pressure distribution. The flow proper- 
ties may also be used to map laminar and turbulent flow 
regions near the surface of the wing and to predict the forma- 35 
tion of shock waves in transonic and supersonic flow. 

A computer-generated simulation can be performed on a 
computer-generated aircraft surface to simulate the fluid 
dynamics of a surrounding fluid flow. The geometry of the 
computer-generated aircraft surface is relatively easy to 40 
change and allows for optimization through design iteration 
or analysis of multiple design alternatives. A computer-gen- 
erated simulation can also be used to study situations that may 
be difficult to reproduce using a physical model, such as 
supersonic flight conditions. A computer-generated Simula- 45 
tion also allows a designer to measure or predict fluid-flow 
properties at virtually any point in the model by direct query, 
without the difficulties associated with physical instrumenta- 
tion or data acquisition techniques. 

In some cases, a computer-generated simulation includes a 50 
computational fluid dynamics (CFD) simulation module used 
to predict the properties of the fluid flow. A CFD simulation 
module estimates the properties of a simulated fluid flow by 
applying a field equation that estimates the interaction 
between small simulated fluid volumes, also referred to as 55 
fluid cells. Because a single CFD simulation module may 
include millions of individual fluid cells, the complexity of 
the relationship between fluid cells can have a large effect on 
the computational efficiency of the simulation. Complex 
CFD simulation modules can be computationally expensive 60 
and require hours or even days to execute using high-perfor- 
mance computer processing hardware. 

To reduce the computational burden, in some instances it is 
desirable to use a CFD simulation module that simplifies the 
fluid dynamics and produces a fluid simulation that can be 65 
solved more rapidly. For example, for fluid flows that are 
relatively uniform or are located away from an aircraft sur- 


2 

face, a simplified simulation that minimizes or ignores certain 
fluid dynamic phenomena can be used. In the examples dis- 
cussed below, a simplified simulation may ignore fluid 
dynamic contributions due to fluid viscosity, which, in some 
cases, have little effect on the overall behavior of the fluid 
flow. A simplified simulation that ignores fluid viscosity may 
be called an inviscid simulation. By using an inviscid simu- 
lation to simulate at least part of the fluid flow, processing 
time may be improved. 

In other situations, where the fluid flow is not as uniform, it 
may be necessary to use a CFD simulation module that is 
more sophisticated and capable of accurately predicting the 
fluid properties, using more complex fluid dynamics. In the 
examples discussed below, a more sophisticated simulation 
may account for dynamic contributions due to fluid viscosity. 
A simulation that accounts for fluid viscosity may be called a 
viscous simulation. Under certain conditions, a viscous simu- 
lation may be more accurate, particularly for portions of the 
fluid flow near the aircraft surface where fluid viscosity 
affects the results. However, viscous simulations are also 
likely to require more computing resources and therefore 
require more time to solve. 

It may be advantageous to construct a hybrid computer- 
generated simulation that employs both an inviscid CFD 
simulation module in locations where the fluid flow is rela- 
tively uniform, and a viscous CFD simulation module in 
locations where the fluid dynamics are more complex. By 
combining different CFD simulation modules, a hybrid com- 
puter-generated simulation may increase processing speed 
while producing accurate results. 

Using multiple CFD simulation modules may be difficult, 
particularly if the CFD simulation modules were not config- 
ured to work together. The interface between the simulation 
modules must be constructed so that the resulting computer- 
generated simulation is both computationally efficient and 
analytically robust. The techniques described herein solve 
some of the difficulties in implementing a computer-gener- 
ated simulation using multiple simulation modules. 

SUMMARY 

One exemplary embodiment includes a computer-imple- 
mented method of generating a fluid-flow simulation over a 
computer-generated surface using a quasi-simultaneous tech- 
nique. The simulation includes an inviscid fluid-flow mesh 
comprised of a plurality of inviscid fluid cells and a viscous 
fluid-flow mesh comprised of a plurality of boundary-layer 
fluid cells. At least some of the boundary-layer fluid cells are 
on or adjacent to the computer-generated surface. 

An initial fluid property is determined for at least one 
inviscid fluid cell using an inviscid fluid simulation that does 
not simulate fluid viscous effects. An initial boundary-layer 
fluid property is determined for at least one of the boundary- 
layer fluid cells using the initial fluid property and a viscous 
fluid simulation that simulates fluid viscous effects. The at 
least one inviscid fluid cell is close in proximity to the at least 
one boundary-layer fluid cell. 

An updated boundary-layer fluid property is determined 
for the at least one boundary-layer fluid cell using the initial 
fluid property, initial boundary-layer fluid property, and an 
interaction law. The interaction law approximates the inviscid 
fluid simulation using a matrix of aerodynamic influence 
coefficients computed using a two-dimensional surface panel 
technique and a fluid-property vector. In some embodiments, 
the interaction law approximates the inviscid fluid simulation 
using a matrix of aerodynamic influence coefficients com- 
puted using an axisymmetric surface panel technique instead 
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of a two-dimensional surface panel technique. An updated 
fluid property for the at least one inviscid fluid cell is deter- 
mined using the updated boundary-layer fluid property. 

In some embodiments, the matrix of aerodynamic influ- 
ence coefficients is adapted to account for compressibility by 5 
dividing each row by (1-Mj 2 ) 172 , wherein Mj is the local 
Mach number. In some embodiments, the matrix of aerody- 
namic influence coefficients is adapted row -by-row to imple- 
ment either the subsonic or supersonic interaction law 
depending on the local Mach number. The matrix of aerody- 10 
namic influence coefficients may be further adapted to utilize 
a linear weighted average of the subsonic and supersonic 
interaction laws to approximate the changes to the inviscid 
fluid flow for local Mach numbers near Mach 1 . 

15 

DESCRIPTION OF THE FIGURES 

FIG. 1 depicts a computer-generated fluid flow applied to a 
computer-generated aircraft surface. 

FIG. 2 depicts an exemplary process for coupling an invis- 20 
cid and boundary-layer CFD simulation module using a 
quasi-simultaneous technique. 

FIGS. 3 a and 3b depict exemplary data flows between an 
inviscid and boundary-layer CFD simulation module. 

FIG. 4 depicts an exemplary surface mesh and a fluid-flow 25 
mesh. 

FIG. 5 depicts an exemplary surface mesh and a fluid-flow 
mesh. 

FIGS. 6 a and 6b depict an exemplary fluid flow around a 
wing surface. 30 

FIG. 7 depicts an exemplary computer system for simulat- 
ing a fluid flow over an aircraft surface. 

FIG. 8 depicts an exemplary computer network. 

The figures depict one embodiment of the present inven- 
tion for purposes of illustration only. One skilled in the art will 35 
readily recognize from the following discussion that alterna- 
tive embodiments of the structures and methods illustrated 
herein can be employed without departing from the principles 
of the invention described herein. 

40 

DETAILED DESCRIPTION 

As discussed above, a computer-generated simulation can 
be used to analyze the aerodynamic performance of a pro- 
posed aircraft surface, such as a fuselage, wing, or control 45 
surface. Using known geometry modeling techniques, a com- 
puter-generated aircraft surface that represents the outside 
surface of the proposed aircraft can be constructed. FIG. 1 
depicts an exemplary computer-generated aircraft surface of 
the Space Shuttle orbit er vehicle, external tank, and twin solid 50 
rocket boosters. A CFD fluid simulation module has been 
applied using the computer-generated aircraft surface of the 
Space Shuttle orbiter to predict the fluid properties of an 
exemplary fluid flow. 

As shown in FIG. 1, the results of the simulation can be 55 
visually represented as shaded regions on the computer-gen- 
erated aircraft surface of the Space Shuttle. Different shades 
represent the predicted pressure distribution resulting from 
the simulated fluid flow. In FIG. 1, transitions between the 
shaded regions represent locations of predicted pressure 60 
change across the surface of the Space Shuttle. Similarly, 
different pressures in the surrounding fluid flow are repre- 
sented as differently shaded regions. 

In FIG. 1, the simulation of the fluid flow is visualized by 
depicting the predicted surface-pressure distribution. How- 65 
ever, the simulation may be visualized using other fluid prop- 
erties, including velocity, air temperature, air density, and 
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others. Additionally, the simulation may be used to visualize 
locations of developing shock waves or transitions between 
laminar and turbulent flow. 

The simulation allows the designer or engineer to evaluate 
the performance of the aircraft geometry for various flow 
conditions. If necessary, changes can be made to the aircraft 
geometry to optimize performance or eliminate an unwanted 
aerodynamic characteristic. Another simulation can be per- 
formed using the modified geometry and the results can be 
compared. To allow for multiple design iterations, it is advan- 
tageous to perform multiple simulations in a short amount of 
time. However, as described above, there is a tradeoff 
between speed and accuracy of the simulation depending on 
the type of CFD simulation module used. 

As described in more detail below, a computer-generated 
simulation can be constructed that uses multiple CFD simu- 
lation modules to simplify the calculations required to 
achieve a complete simulation of the fluid flow. The calcula- 
tions can be further simplified by using an interaction law to 
approximate the results of one of the CFD simulation mod- 
ules when perfonning the simulation. Specifically, an inter- 
action law allows a first, viscous CFD simulation module to 
determine one or more predicted fluid properties by using the 
interaction law to estimate the predicted fluid properties of a 
second, inviscid CFD simulation module. For CFD simula- 
tion modules with results that are highly interdependent, an 
interaction law can be used to improve the speed and reliabil- 
ity of the simulation results. 

As mentioned above, a computer-generated simulation 
may represent a fluid flow using a mesh of fluid cells. The size 
and shape of the fluid cells can vary depending on the type of 
fluid-flow mesh used. In some cases, the fluid cells are small, 
cuboid volumes of fluid. In other cases, the fluid cells are 
small, thin volumes of fluid, each thin volume having only a 
nominal thickness. In other cases, the fluid cells may be 
represented by a two-dimensional polygon or a single point. 

In the examples provided below, inviscid regions of the 
fluid flow (typically away from the aircraft surface) are rep- 
resented using a fluid-flow mesh of three-dimensional, 
cuboid fluid cells. Viscous regions of the fluid flow (typically 
near the surface of the aircraft) are represented using a fluid- 
flow mesh of either three-dimensional fluid cells, two-dimen- 
sional fluid cells, or a set of boundary-layer prediction points 
located along the computer-generated aircraft surface. 

Both an inviscid and viscous CFD simulation modules 
predict the interactions between the fluid cells in the fluid- 
flow mesh, using a fundamental algorithm or field equation. 
The speed and accuracy of a CFD simulation module 
depends, in part, on the field equation used to predict the 
interaction between the flow cells. For an inviscid simulation, 
the field equation simplifies the relationship between flow 
cells by ignoring or minimizing viscous contributions. Vis- 
cous effects of a fluid can be minimized or ignored when, for 
example, there is not a significant velocity difference between 
adjacent fluid cells, and therefore shear forces due to internal 
friction or viscosity are minimal compared to inertial forces. 
Field equations for an inviscid simulation are typically less 
complex and require a less dense fluid-flow mesh, and there- 
fore are more computationally efficient. F or instance, an algo- 
rithm that solves the Euler equations may be used to simulate 
an inviscid fluid flow. 

In other instances, a more complex field equation is used to 
more accurately predict the interaction between the fluid 
cells. For example, a method that solves the Navier-Stokes 
equations can be used to simulate the pressure and shear 
forces on the fluid cells. Unlike the Euler solver mentioned 
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above, a Navier-Stokes solver accounts for the effects of 
viscosity and offers a more accurate simulation of a fluid flow. 

However, the improved accuracy of a Navier-Stokes-based 
solution comes at the cost of increased computational load 
and increased operator preparation time, and therefore a 5 
Navier-Stokes solver is generally slower than an Euler-based 
algorithm. Thus, selecting the field equation for a CFD mod- 
ule often involves a tradeoff between speed and accuracy. In 
practice, designers may use faster Euler-based CFD models 
to evaluate multiple design iterations and then validate the to 
final design iteration with a more accurate Navier-Stokes- 
based CFD model. However, if the Navier-Stokes CFD simu- 
lation reveals a design problem, the entire process must be 
repeated, wasting valuable time and computing resources. 

The techniques described below are computer-generated 15 
simulations that use both inviscid and viscous field equations 
to achieve acceptable accuracy without requiring the compu- 
tational burden of a full viscous CFD simulation (e.g., a 
Navier-Stokes-based CFD simulation). In many simulations, 
there is a region of the fluid flow that can be accurately 20 
predicted without taking viscous contributions into account. 

For example, regions of the fluid flow that are located away 
from an aircraft surface or wing surface, have a relatively 
uniform velocity profile. Therefore, an inviscid simulation 
using, for example, an Euler-based analysis, can be used to 25 
accurately predict the behavior of the fluid-flow region. In 
other fluid-flow regions, where there is a less uniform velocity 
profile, a more complex, viscous simulation can be used. 

It may be advantageous to generate a simulation using 
more than one field equation to simulate the fluid flow over a 30 
computer-generated aircraft surface. In particular, it may be 
advantageous to generate a simulation using an inviscid CFD 
simulation module (using an inviscid field equation) to simu- 
late regions of the flow away from the aircraft surface and a 
viscous CFD simulation module (using a viscous field equa- 35 
tion) to simulate regions of the flow near the aircraft surface. 
This technique reduces the computational load by simplify- 
ing the simulation for regions of the flow that do not need the 
more complex, viscous CFD simulation. 

Typically, the inviscid (intertia-dominated) and viscous 40 
(viscosity-dominated) regions of the fluid flow are in physical 
contact with each other, and therefore the predicted fluid-flow 
properties of, for example, the inviscid region may influence 
the fluid properties of the viscous region. In some cases, there 
is a strong interaction between the inviscid and viscous 45 
regions of the fluid flow, and small changes in the inviscid 
fluid properties may result in a large change in the viscous 
fluid properties, and vice versa. This is particularly true near 
trailing edges of wings and where the fluid-flow is transition- 
ing from supersonic to subsonic flow. 50 

A simulation that uses multiple CFD simulation modules 
to simulate the different regions of the fluid flow must manage 
the interface between the CFD simulation modules to accu- 
rately simulate the interdependencies between the respective 
simulated regions of the fluid flow. Typically, this requires 55 
that the multiple CFD modules converge to the same result 
even though the CFD modules use different field equations to 
predict the fluid properties. 

In one example, fluid properties predicted for each region 
of the fluid flow are passed between the respective CFD 60 
simulation modules. Sometimes referred to as the direct 
method, the fluid properties of, for example, an inviscid CFD 
simulation module are passed to viscous CFD simulation 
module. In turn, the results of the viscous CFD simulation 
module are passed back to the inviscid CFD simulation mod- 65 
ule. This is easily accomplished because the natural output of 
the inviscid CFD simulation module is the natural input of the 


6 

viscous CFD simulation module, and, in turn, the natural 
output of the viscous CFD simulation module is the natural 
input of the inviscid CFD simulation module. Examples of 
natural inputs and outputs of inviscid and viscous CFD simu- 
lation modules are discussed in more detail below with 
respect to equations 13, 18, 19, and 20. These natural forms of 
the CFD simulation modules are called direct solvers, leading 
to them being designated as a “direct method.” This exchange 
is conducted for multiple iteration steps until a convergent 
solution is determined. It is also possible to create inverse- 
mode versions of the inviscid and viscous CFD solvers where 
the inputs and outputs are swapped. The inverse solvers can 
be used to formulate the direct-inverse (a.k.a. semi-inverse) 
method and the inverse-inverse (a.k.a. full-inverse) method 
for solving the coupled viscous-inviscid problem. 

During each iteration step, the CFD simulation modules 
determine the fluid properties for a respective region of the 
fluid flow without accounting for dynamic contributions from 
other regions of the fluid flow. As a result, there may be 
drawbacks to performing a direct method simulation. In some 
cases, the two CFD simulation modules may not converge to 
a single solution. In other cases, the CFD simulation modules 
may converge but so slowly as to not be of practical utility. 

As an alternative, the field equations from the CFD simu- 
lation modules may be combined to determine a single set of 
fluid-properties for a combined fluid-flow region. That is, 
both the inviscid and viscous field equations are solved simul- 
taneously so that the dynamic contributions from the two 
simulated regions of the fluid flow are accounted for. This 
technique may also be called a simultaneous simulation. 
There are, however, drawbacks to performing a simultaneous 
simulation. First, constructing a simultaneous viscous-invis- 
cid simulation is a very complex task and has rarely been 
accomplished for three-dimensional flow simulations. Sec- 
ond, the simultaneous simulation technique cannot be applied 
using existing inviscid CFD simulation modules. Instead, a 
simultaneous simulation requires specially-crafted inviscid 
and viscous simulation modules that can work with a single 
equation solver capable of simultaneously solving many non- 
linear equations. This impairs the compatibility with existing 
simulation tools, and therefore may dramatically increase 
development and implementation costs. 

As discussed in more detail below, inviscid and viscous 
CFD simulation modules can be coupled using a quasi-simul- 
taneous technique. Using a quasi-simultaneous technique, 
the viscous CFD simulation module uses an interaction law to 
account for dynamic changes to the inviscid fluid properties 
due to changes in the viscous region of the fluid flow. 

Using a quasi-simultaneous technique, an existing inviscid 
CFD simulation module can be used to predict a set of initial 
inviscid fluid property values (for an inviscid region of fluid 
flow). A set of viscous fluid properties (for a viscous region of 
the fluid flow) are determined using the initial inviscid fluid 
properties, a viscous CFD simulation module, and an inter- 
action law that approximates changes to the inviscid fluid 
properties due to presence of the viscous region of the fluid 
flow. 

A quasi-simultaneous simulation may provide several 
advantages. First, because the viscous CFD simulation mod- 
ule includes a prediction of the changing inviscid fluid prop- 
erties, iterated results between the two CFD simulation mod- 
ules may converge more rapidly. Additionally, by selecting an 
interaction law that enforces key boundary conditions, the 
converged results are more likely to result in a robust and 
accurate prediction of the fluid flow. Specifically, an interac- 
tion law based on linearized potential flow equations can be 
used to enforce the Kutta condition at the trailing edge of the 
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wing surface, and therefore ensures that fluid-flow informa- 
tion is propagated between the top and bottom viscous flow 
layers near the computer-generated aircraft surface (particu- 
larly near a computer-generated wing surface). 

The techniques described herein provide a quasi-simulta- 5 
neous simulation that combines an inviscid and a viscous 
CFD simulation module using an interaction law based on 
linearized potential flow. The following discussion provides 
an example of a simulated fluid flow over a wing surface. 
However, the technique may also be applied to a simulated 10 
fluid flow over any type of surface subjected to a fluid flow. 

For example, the following techniques could be applied to the 
surface of a space vehicle, land vehicle, watercraft, or other 
object having a surface exposed to a fluid flow. In addition, the ; 5 
following techniques can be applied to simulations of various 
types of fluid flow, including, for example, a gas fluid flow or 
liquid fluid flow. 

1 . Quasi-Simultaneous Technique for Combining CFD Simu- 
lation Modules 20 

As discussed above, a computer-generated simulation can 
use both inviscid and viscous CFD simulation modules to 
represent different regions of the fluid flow. The two CFD 
simulation modules must be linked so that the results from 
each module provide an accurate and realistic prediction of 25 
the combined fluid-flow regions. 

In general, and as explained in more detail below in Section 
4, an inviscid CFD simulation module can be used to simulate 
outer fluid-flow regions that are at least some distance away 
from an aircraft surface. Also explained in Section 4, a vis- 30 
cous CFD simulation module can be used to simulate fluid- 
flow regions that are adjacent to an aircraft surface, such as a 
boundary-layer fluid-flow region close to a wing surface. 

FIG. 2 depicts an exemplary process 200 for using a quasi- 
simultaneous technique to couple inviscid and viscous CFD 35 
simulation modules. More specifically, the exemplary pro- 
cess 200 uses quasi-simultaneous technique with an interac- 
tion law to approximate outer (inviscid) fluid-flow effects on 
the simulated boundary-layer (viscous) fluid- flow region and 
vice versa. As explained in more detail below, an interaction 40 
law that uses linearized potential flow can be used to enforce 
key boundary conditions on the simulation so that the pre- 
dicted result is both an accurate and realistic prediction of the 
combined fluid flow. 

In the following discussion, for simplicity, the result of the 45 
inviscid simulation is assumed to depend only on the two- 
dimensional displacement thickness, 8*, and return the 
boundary-layer edge velocity, u(u„=E(8*„)). In general, how- 
ever, the result of an inviscid simulation may depend on the 
geometry of the entire aircraft including the surface slope of 50 
both the chordwise and spanwise displacement thicknesses in 
three dimensions (or instead on the transpiration flux which 
itself depends on the two displacement thickness slopes). 
Similarly, the boundary-layer simulation is simplified to 
depend only on the boundary-layer edge velocity u and to 55 
return the displacement thickness (8* ?i =B(u„)). However, in 
general, the boundary-layer simulation may depend on the 
many fluid properties at the boundary-layer edge (velocity 
vector, pressure, density, etc.) and can return any number of 
boundary-layer parameters including the chordwise and 60 
spanwise displacement thicknesses, transpiration flux, and 
transpiration velocity. 

In operation 202, an initial fluid property, for example 
velocity u„, is determined for one or more inviscid fluid cells 
420 of a fluid-flow mesh. As explained in more detail below 65 
with respect to FIGS. 4 and 5, an inviscid fluid flow can be 
simulated using a fluid-flow mesh 402 of inviscid fluid cells. 
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The fluid properties of each inviscid fluid cell (including the 
initial fluid velocity u„) are predicted using an inviscid CFD 
simulation module. 

The initial fluid velocity u„, also called the edge fluid 
velocity, may be also expressed as a solution of the Euler 
equations E(8*„), which may be dependent on the initial 
boundary-layer thickness 8*„. In some cases, the initial fluid 
velocity u„ is a vector quantity represented using three-di- 
mensional velocity components. 

An exemplary inviscid CFD simulation module is dis- 
cussed below with respect to inviscid field equation 13. The 
initial fluid velocity u„ may be determined using equation 13 
and a set of initial conditions used to seed the simulation. The 
inviscid CFD simulation module may iterate multiple times 
before determining the initial fluid velocity value u„. In some 
cases, the initial fluid velocity u„ represents the outer (invis- 
cid) fluid-flow conditions ignoring effects due to a boundary- 
layer fluid-flow region. 

In operation 204, initial boundary-layer fluid properties, 
for example, boundary-layer thickness 8*„ is determined for 
at least one boundary -layer fluid cell 430. The initial bound- 
ary -layer thickness 8*„ is determined using equations 18, 19, 
and 20 and the inviscid fluid properties, including the initial 
fluid velocity u„. The initial boundary-layer thickness 8*„ 
may be a predicted thickness or transpiration flux, as dis- 
cussed in more detail below with respect to equations 18, 19, 
and 20. Hie initial boundary -layer fluid velocity u w „ can also 
be determined using boundary -layer field equations 18 and 

1 9, if solved in inverse mode. The initial boundary-layer fluid 
velocity u w n can also be expressed as B(8*„), which is depen- 
dent on initial boundary -layer thickness 5*„. 

As shown in FIGS. 4 and 5, the at least one boundary -layer 
fluid cell 430 should be close in proximity to one or more 
inviscid fluid cells 420. The fluid cells 430 and 420 should be 
proximate so that the updated boundary-layer thickness 8*„ +1 
can be influenced by of the nearby regions of the inviscid fluid 
flow. 

In operation 206, an updated boundary-layer fluid property 
8*„ +1 is determined for the at least one boundary-layer fluid 
cell using the boundary-layer equations and an interaction 
law 1(8*) to approximate the inviscid fluid simulation. The 
updated boundary-layer fluid property 8*„ +1 may be sensitive 
to changes in the inviscid fluid-flow region, particularly if 
there is a strong interaction between the inviscid and bound- 
ary-layer fluid-flow regions. An interaction law 1(8*) 
approximates changes in the inviscid fluid flow due to 
changes in the boundary-layer fluid flow. A more detailed 
discussion of an exemplary interaction law is included in 
Section 2, below. 

Using a quasi-simultaneous technique, an updated bound- 
ary-layer fluid property value S* M+1 can be determined by 
simultaneously solving for the boundary-layer fluid proper- 
ties at the time n+ 1 (using, for example, equations 18,19, and 

20, below) and a prediction of the inviscid fluid-flow proper- 
ties at time n+1 (using, for example, 1(8*) in equation 1, 
below). 

In general, the prediction of an inviscid fluid velocity at 
time n+1 (u'„ +1 ), can be expressed in terms of an initial invis- 
cid fluid velocity at time n(u„ or E(8*„)) and an interaction 
law 1(8*): 

»'„ tl =£(6*„)+/(6*„ tl )-/(6*„) Equation 1 

As discussed in more detail below with respect to equations 5 
and 7 below, using an interaction law 1(8*) that is based on the 
principles of linear potential flow, the results may more accu- 
rately reflect real fluid dynamics by enforcing boundary con- 
ditions, including the Kutta condition. Additionally, because 
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the interaction law 1(8*) is only used to update the inviscid 
fluid velocity due to an incremental change in the displace- 
ment thickness since the previous iteration, it does not affect 
the converged solution, and therefore, the absolute accuracy 
of interaction law 1(8*) prediction is not critical. 5 

Using a quasi-simultaneous coupling technique, the inter- 
action law 1(8*) can be used as a proxy for updated fluid 
properties that would otherwise be determined using the 
inviscid CFD simulation module. As shown in equation 1 , the 
interaction law 1(8*) can be used to predict the inviscid fluid 10 
velocity at time n+1 (u'„ +1 ) to estimate changes to the inviscid 
fluid-flow region without having to perform another complete 
simulation using the inviscid CFD simulation module. The 
relative simplicity of the interaction law 1(8*), allows the 
interaction law to be combined with boundary-layer field 
equations 1 8 and 1 9 to obtain an improved prediction of the 
boundary-layer fluid properties. 

As discussed above, a boundary-layer CFD simulation 
module (using boundary-layer field equations 18 and 19, 
below) can be used to determine an initial boundary-layer 
fluid velocity u„ bl given an initial boundary-layer displace- 
mentthickness 8*„ (thus, u„ W =B(8*„)). Similarly, an inviscid 
CFD simulation module can compute a fluid velocity u„ given 
an initial boundary-layer displacement thickness 8*„ (thus, 
u„=E(8*„)). Specifying that u=u w and using Newton’s 
method to solve the simultaneous equations for an updated 
boundary-layer thickness 8*„ +1 , a Newton iteration can be 
expressed as: 


dS * dS- 


Equation 2 


35 

Similarly if equation 1 is used to solve for u' H+1 simulta- 
neously with u'„ +1 =u„ +1 =B(8*„ +1 ) using Newton’s method 
the updated boundary-layer thickness 8*„ +1 can be expressed 
as: 


dl 

IF 


dB 

' IF 


(B(S* n )-E(SD). 


Equation 3 


45 

Expressed in terms of fluid velocities u bln and u,„ the updated 
boundary-layer thickness 8*, J+1 can be written as: 


1 Equation 4 

s ^ =s " + -dl— 1FFF: (Ubln ^ Unl 

dS * dS- 

In this way, an updated boundary-layer thickness S* ;;+1 can 55 
be determined using initial boundary-layer thickness 8*„, the 
initial inviscid fluid velocity u„, initial boundary-layer fluid 
velocity u„ b! , and gradient of the interaction law 


dS 

The advantage is that the gradient of the full inviscid equa- 
tions (dE/d8*) is generally too complicated to be calculated in 65 
a practical manner (as would be needed in order to apply the 
simultaneous method mentioned above). Equations 2 through 


4 illustrate that the quasi-simultaneous method lias some 
qualities in common with the simultaneous method solved 
using an approximate Newton method. Thus the designation 
as a “quasi-simultaneous” method. 

In operation 208, an updated fluid velocity u„ +1 orE(8*„ +1 ) 
is determined for the one or more inviscid fluid cells using the 
updated boundary-layer thickness 8*„ +1 . The updated bound- 
ary-layer thickness 8*„ +1 may affect the inviscid or outer fluid 
flow by increasing or decreasing the size of the boundary- 
layer region effectively changing the shape of the inviscid 
fluid-flow region around the aircraft surface. An updated fluid 
velocity u„ +1 can be calculated using an inviscid CFD simu- 
lation module (using, for example, equation 13, below). The 
updated fluid velocity u H+1 or E(8* n+1 ) for the one or more 
inviscid fluid cells can be used to update the other inviscid 
fluid cells in the fluid-flow mesh. 

The operations of process 200 can be repeated using the 
updated inviscid fluid u„ +1 as a new initial inviscid fluid veloc- 
ity u„ (above). In this way, process 200 can be iterated mul- 
tiple times until convergent results are achieved. 

FIGS. 3 a and 3 b depict an exemplary data flow 300 for 
iterating results between an inviscid CFD simulation module 
310 and a viscous or boundary-layer CFD simulation module 
320. The data flow 300 can be used to link an existing inviscid 
CFD simulation module 320 with a viscous CFD simulation 
module 310 that has been adapted to use a quasi-simultaneous 
technique. 

In operation 312, an initial fluid velocity u ; (u„ where n=l) 
is determined for one or more inviscid fluid cells 420 of a 
fluid-flow mesh. As explained above with respect to operation 
202, the fluid properties of each inviscid fluid cell (including 
the initial fluid velocity iq) can be predicted using the inviscid 
CFD simulation module 310. Specifically, the initial fluid 
velocity iq may be determined using equation 13 and a set of 
initial conditions used to seed the simulation. The initial fluid 
velocity iq is passed to a boundary-layer CFD simulation 
module 320. 

In operation 322, initial boundary-layer fluid properties, 
boundary-layer thickness b* 1 or velocity u w x (8*„ and u w „ 
where n=l) are determined for at least one boundary-layer 
fluid cell 430. As explained above with respect to operation 
204, the initial boundary -layer thickness 8* : or initial bound- 
ary-layer fluid velocity u w j can be determined using bound- 
ary-layer field equations 18, 19, and 20 of the boundary-layer 
CFD simulation module 320. 

In operation 324, a change to the inviscid fluid properties is 
estimated by applying equation 1 . Discussed in more detail in 
Section 2 below, an interaction law that uses linearized poten- 
tial flow can be used to enforce key boundary conditions on 
the simulation. 

In operation 326, an updated boundary-layer fluid property 
8*„ +1 is determined for the at least one boundary-layer fluid 
cell using the result of operation 324, which uses an interac- 
tion law 1(8*) to approximate the inviscid fluid simulation, as 
explained above with respect to operation 206. The updated 
boundary-layer fluid property 8* n+1 is passed to the inviscid 
CFD simulation module 310. 

As shown in FIGS. 3a and 3 b, the operations 324 and 326 
may be repeated until the iteration converges to a compatible 
solution, or, preferably, operations 324 and 326 are solved 
simultaneously using standard iterative techniques for solv- 
ing multiple nonlinear equations such as Newton’ s method or 
Gauss-Seidel relaxation. 

In operation 316 an updated fluid velocity u„ +1 or E(8*„ +1 ) 
is determined for the one or more inviscid fluid cells. As 
explained above with respect to operation 208, the updated 
fluid velocity u„ +1 can calculated using an inviscid CFD simu- 
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lation module 310 (using, for example, equation 13, below). 
The updated fluid velocity u M+1 or E( 8 *„ +1 ) for the one or 
more inviscid fluid cells can be used to update the other 
inviscid fluid cells in the fluid-flow mesh. 

As shown in FIG. 3a, the updated fluid velocity u „ +1 can 
also be used to seed another iteration of data flow 300 with the 
updated fluid velocity u M+1 serving as the fluid velocity u„ in 
operation 326. Alternatively, as shown in FIG. 3 b. the updated 
fluid velocity u M+1 can be used to seed another iteration of data 
flow 300 with the updated fluid velocity u „ +1 serving as the 
fluid velocity u„ in operation 322. In this way, data flow 300 
can be iterated multiple times until convergent results are 
achieved. 

2. Exemplary Interaction Law 

As discussed above, an interaction law can be used to 
estimate updated inviscid fluid properties due to changes in a 
nearby viscous fluid flow. A viscous CFD simulation module 
incorporating an interaction law can anticipate changes to the 
inviscid fluid flow without having to perform a complete 
inviscid simulation using, for example, an inviscid CFD 
simulation module. For neighboring fluid flows that are 
highly interdependent, an interaction law can be used to 
improve the speed and reliability of the simulation results. 

As mentioned above, an interaction law can also be used to 
enforce boundary conditions so that the physical constraints 
of an actual fluid flow over and around a wing surface are 
maintained. For example, the Kutta condition at the trailing 
edge of an aircraft surface is enforced ensuring that the fluid- 
flow information is propagated between, for example, the top 
and bottom surfaces of a computer-generated wing surface. 

In the example below, a two-dimensional surface panel 
technique using linearized potential flow theory is used as the 
basis for the interaction law. This approach simplifies the 
calculation by representing the inviscid fluid flow as a two- 
dimensional strip along a cross section of a computer-gener- 
ated aircraft surface. The following example is performed for 
a computer-generated wing surface (as shown in FIGS. 6 a 
and 6b). 

For a computer-generated wing surface the interaction law- 
can be expressed as: 

lib* )=ld.jmj. Equation 5 

where the variable m 7 represents the product of boundary- 
layer fluid properties: edge density p c , edge velocity u e , and 
boundary-layer thickness 8 *: 

m=p L> u t> 6*. Equation 6 

Similarly an interaction law for a swept/tapered wing sur- 
face can be expressed as: 

7(6* )=lb.jV^. Equation 7 

where v w is the transpiration velocity at the surface. Matrices 
d, ; and h tJ are fimctions of the computer-generated aircraft 
surface geometry and are also referred to as aerodynamic 
influence coefficient matrices. 

The interaction influence coefficients in equations 5 and 7 
can be computed using a two-dimensional surface panel tech- 
nique. For example, matrices d ( / and b l f can be derived using 
a set of flat panels containing sources of constant strength and 
vortices of linearly varying strength that shape the approxi- 
mated inviscid fluid flow. The sources induce a transpiration 
flow perpendicular to the computer-generated aircraft sur- 
face, thus modeling the effect of a boundary-layer fluid-flow 
region having a thickness 8 *. Collectively, the sources and 
vortices modify the path of an approximated inviscid fluid- 
flow region to flow around the computer-generated aircraft 
surface and boundary-layer fluid-flow region (having a thick- 


ness 8 *). An inviscid fluid velocity u can be estimated using 
this approximated inviscid fluid flow. As suggested above, 
this approach can also be used to enforce the Kutta condition 
at the trailing edge of a computer-generated wing surface, and 
5 therefore ensures that flow information is propagated from 
the region of the flow near the trailing edge of the wing 
surface to the flow around other portions of the wing surface 
encompassing both the top and the bottom surface of the 
computer-generated wing geometry. Moreover, specific coef- 
111 ficient schemes can be implemented in matrices d, and b, to 
account for either blunt or sharp trailing edges. 

Equations 5 and 7 apply only for subsonic incompressible 
flow around a wing surface when the influence coefficients 
15 are generated using a linear surface panel technique, as 
described above. To increase applicability to higher-speed 
flows up to transonic flows, a compressibility correction can 
be added to the interaction law following the transonic small 
disturbance corrections described by Spreiter and Alksne. 
20 See, e.g., Spreiter, J. R. and Alksne, A. Y., “Thin Airfoil 
Theory Based on Approximate Solution of the Transonic 
Flow Equation,” NASA TR-1359, 1958. Using this tech- 
nique, the magnitude of each row of the b tJ or A tJ matrix can 
be increased according to the local sweep-perpendicular 
25 Mach number by dividing by the tenn (1 -Mj 2 ) 1 ' 2 . The value 
of the Mach number is derived from the solution of the Euler 
equations (for example, equation 13). Here, sweep-perpen- 
dicular Mach number is the component of local velocity in a 
particular direction divided by the local speed of sound. That 
30 particular direction may be taken to be: perpendicular to the 
local geometric wing sweep; along a set of points, volumes, or 
fluid cells on which the boundary-layer equations are solved; 
or perpendicular to the local sweep of the lines of constant 
pressure (isobars) on the computer-generated surface. 

35 Other compressibility corrections may be used. In two- 
dimensional flows, the so-called Karman-Tsien correction 
may be applied. In three-dimensional flows, the Prandtl- 
Glauert correction may be applied, although the Prandtl- 
Glauert correction may be less accurate. The improved accu- 
40 racy of the Karman-T sien over Prandtl-Glauert may be due, in 
part, to a better approximation of the local Mach number. For 
fluid simulations where the local Mach number cannot be 
calculated directly, an approximation may be made. How- 
ever, when the outer (inviscid) fluid flow is computed using 
45 the Euler equations (e.g., equation 13) the local Machnumber 
is available directly. 

For supersonic flow condition, equations 8 and 9, below, 
can be used to compute edge velocity u e in the interaction law. 
Equation 8 can be used for a two-dimensional strip boundary- 
50 layer estimate: 


1 ( dni( n ) dm( n+ i ) Equation 8 

55 

and equation 9 can be used for three-dimensional and for 
sweep/taper-based boundary-layer estimates: 

60 

1 Equation 9 

u r<>. I I ) = (V*tn+ 1) “ 

V - 1 

65 In some cases the matrices d,- . and b,- ■ can be adapted row by 

row to implement either the subsonic or supersonic interac- 
tion law depending on the local Mach number. Near Mach 



US 8,892,408 B2 


13 

1.0, the interaction laws may not work. The fluid flow is 
extremely nonlinear and a simple interaction law cannot pre- 
dict the flow behavior. A reasonable work-around is to instead 
utilize a linear weighted average of the subsonic and super- 
sonic interaction laws to approximate the changes to the 5 
inviscid fluid flow very close to Mach 1 .0. For example, the 
following formula can be used to create a compressible influ- 
ence coefficient that works from subsonic to supersonic 
flows, with linear weighted averaging as described below: 


14 

curvature correction can be used to approximate this effect 
and to improve the accuracy of the interaction law. If the 
streamwise geometric curvature of the aircraft surface is k M „ 
then the curvature k of a fluid-flow streamline at the bound- 
ary-layer edge may be approximated by: 

.(MUM Equation 1 1 


: (1 ~Wi) 


- b;; — W; 


6 ;; 


Equation 10 


l -M} yj Mf - t 


15 


The curvature k induces a difference between the velocity at 
the surface u M , (computed using the inviscid CFD simulation 
module) and the velocity at the boundary-layer edge u e which 
may be approximated by: 


where b comp is the compressible influence coefficient matrix 
that can be used directly in equation 7 in the place olb, 7 . The 
compressible influence coefficient matrix b comp tJ can allow 
equation 7 to be applicable for subsonic, transonic and super- 
sonic flows, b, is the incompressible influence coefficient 
matrix derived from the linear surface panel technique, as 
described above. Here, b tJ represents the Kroneker delta func- 
tion (not to be confused with the displacement thickness 8*) 
which is simply notation describing that the right-hand term 
in the above equation only applies to the diagonal elements of 
the matrix. M, is again the sweep-normal local Mach number 
at the streamwise station i, except that in the left-hand term of 
the above equation, M, is capped above at the value of 0.97 
and the value of M, in the right-hand term is capped below at 
the value of 1.03. The value of w ; . is zero when the actual M, 
is less than 0.97, w, is one when the actual M, is greater than 
1.03, and w, varies linearly between zero and one when the 
actual M, lies between 0.97 and 1 .03 . An analogous technique 
can be performed with A tj to make equation 5 applicable from 
subsonic through supersonic flows. 

The interaction laws discussed above are particularly well 
suited for approximating inviscid fluid-flow conditions 
around a computer-generated wing surface. Other aircraft 
surfaces, such as fuselage-like bodies, including nacelles and 
fuel tanks may require different interaction laws. For 
example, the interaction law for subsonic flow around a fiise- 
lage-like body may consist of an axisymmetric panel code. 

The axisymmetric panel code, which approximates the 
fuselage-like body as a locally axisymmetric body, may use a 
ring-source of constant strength, such as that described by 
Hess and Smith and Lewis. See Hess, J. L. and Smith, A. M. 
O., “Calculation of potential flow about arbitrary bodies,” 
Progress in Aerospace Sciences, pp. 1-138, 1967 and Lewis, 
R. I., “Vortex element methods for fluid dynamic analysis of 
engineering systems,” Cambridge University Press, 2005. 
This implementation differs from the two-dimensional inter- 
action law in that there are no vortex elements. Both the shape 
of the body and the extra displacement thickness due to the 
boundary layer are accounted for by the ring-source elements. 

Additionally, the interaction laws discussed above can be 
further adapted to account for other physical phenomenon, 
such as the so-called curvature correction. Typically, a bound- 
ary layer is assumed to be very thin and the outer inviscid fluid 
solution (determined using, for example, equation 13) can be 
approximated as constant between the actual edge of the 
boundary layer and the surface of the aircraft geometry. How- 
ever, in regions of strong streamwise curvature of the bound- 
ary-layer edge, when, at the same time, the boundary layer is 
relatively thick (such as near wing trailing edges), the inviscid 
fluid solution may vary by a discernable amount between the 
boundary-layer edge and the aircraft surface. In this case, a 


u e -u w =ku e (Q+&*), Equation 12 

where 0 is the momentum thickness and 8* the displacement 
thickness. The relation in equation 1 2 may be added to equa- 
tion 1 to improve the accuracy of the method. 

3. Inviscid and Viscous Fluid-Flow Meshing Techniques 

As discussed above, a computer-generated simulation may 
represent a fluid flow using a mesh of fluid cells. The size and 
_, 5 shape of the fluid cells can vary depending on the type of 
fluid-flow mesh used. In the examples provided below, invis- 
cid regions of the fluid flow (typically away from the aircraft 
surface) are represented using a fluid-flow mesh of three- 
dimensional, cuboid fluid cells. Viscous regions of the fluid 
30 flow (typically near the surface of the aircraft) are represented 
using a fluid-flow mesh of boundary-layer fluid cells near the 
computer-generated aircraft surface. In some cases, the 
boundary-layer fluid-flow mesh is a two-dimensional mesh of 
fluid cells, each fluid cell comprising a single point or two- 
35 dimensional area. 

FIG. 4 depicts a cross-sectional view of an inviscid fluid- 
flow mesh 402 of inviscid fluid cells and a viscous fluid-flow 
mesh 404 of boundary-layer fluid cells, both fluid-flow 
meshes represented as structured meshes. While the cross- 
40 sectional view of FIG. 4 depicts the fluid cells as areas, the 
fluid cells in this example are actually three-dimensional 
volumes. 

FIG. 4 also depicts a surface mesh 406 of quadrilateral 
polygons representing the surface of a wing. The surface 
45 mesh 406 should approximate the curvilinear shape of the 
wing surface without creating gaps or breaks between quad- 
rilateral polygons. In some cases, the surface mesh 406 of 
polygons can be used to define a set of fluid cells, each fluid 
cell comprising a single point or two-dimensional area. For 
50 example, the centroid of one or more of surface mesh poly- 
gons can be used to define one or more fluid cells. 

The inviscid fluid-flow mesh 402 and viscous fluid-flow 
mesh 404 shown in FIG. 4 are shown as structured fluid-flow 
meshes defined using a set of vertices of the surface mesh 406 . 
55 For a set of four adjacent vertices on the surface mesh 406 , a 
volume 410 is projected from the surface of the wing in a 
direction as close to a surface normal as possible. The volume 
is partitioned into fluid cells 408 by defining at least two 
surfaces 412 that offset a given distance from the surface of 
60 the wing. 

FIG. 5 depicts another technique for creating surface and 
fluid-flow meshes. In FIG. 5 , the inviscid fluid-flow mesh 502 
is depicted as being a Cartesian mesh. A Cartesian mesh is 
defined as a mesh of cube or rectangular cuboid fluid cells 
65 508 . That is, each fluid cell 508 is bounded by six flat faces 
where opposite faces are parallel and adjacent faces are 
orthogonal. In some cases, laiger fluid cells are divided into 
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smaller fluid cells by defining additional faces at the midpoint 
of the existing faces. Thus, the fluid cells in a Cartesian mesh 
can be different sizes. 

FIG. 5 depicts one example of a boundary-layer fluid-flow 
mesh 504 . This fluid-flow mesh can be created using a struc- 5 
tured mesh as described above with respect to FIG. 4 . Alter- 
natively, the boundary-layer fluid-flow mesh 504 can be con- 
structed as a Cartesian mesh similar to the inviscid fluid- flow 
mesh 502 . 

FIG. 5 also depicts a surface mesh 506 that is constructed 10 
using triangular mesh elements rather than quadrilateral mesh 
elements. In some cases, the surface mesh 506 of polygons 
can be used to define a set of fluid cells, each fluid cell 
comprising a single point or two-dimensional area. For L - 
example, the centroid of one or more surface mesh polygons 
can be used to define one or more fluid cells. 

4. Field Equations for Simulating Fluid Flow Over a Com- 
puter-Generated Wing Surface 

As discussed above, a fluid flow can be simulated using 20 
both an inviscid and viscous CFD simulation modules to 
improve computing efficiency. The explanation below 
describes how the CFD simulation modules can be applied to 
different regions of fluid flow and provides exemplary field 
equations for simulating the dynamics of the fluid flow. 25 

FIGS. 6 a and 6b depict a cross-sectional representation of 
a fluid flow over a wing surface 602 classified by two regions: 
an outer region 604 and a boundary-layer region 606 . Note 
that the boundary-layer region is drawn, for illustrative pur- 
poses, much thicker than a typical boundary-layer region on 30 
an aircraft wing. Also note that for internal flows, such as for 
flow through a pipe or aircraft engine inlet, the outer region is 
taken to be the region outside the boundary-layer region even 
though that region may be surrounded by geometric surfaces. 

As shown in FIG. 6 a, the outer region 604 may be located 35 
away from the wing surface 602 . However, the outer region 
604 may be close to the wing surface 602 in areas where the 
boundary-layer region 606 is thin or has yet to develop. See, 
for example, the portion of the fluid flow in FIG. 6 a near the 
leading edge of the wing surface 602 . The outer region 604 is 40 
usually characterized as having a relatively uniform velocity 
profile 612 . When there is a relatively uniform velocity profile 
612 , internal shear forces acting on a fluid may be relatively 
small compared to internal forces, and therefore viscous con- 
tributions to the fluid dynamics can be minimized or ignored. 45 
Neglecting viscous effects in a fluid flow results in what is 
known as an inviscid flow. 

A CFD simulation module that ignores viscous effects may 
also be called an inviscid CFD simulation module. Equation 
13, below, provides an exemplary field equation for an invis- 50 
cid CFD simulation module. Equation 13, also called the 
Euler’s equation, represents the conservation of mass, con- 
servation of three components of momentum, and conserva- 
tion of energy: 
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where, u, v, and w are components of the velocity vector, p is 
the pressure, p is the density, and E is the total energy -per-unit 
volume. Combining equation 13 with an equation of state 
(e.g., the ideal gas law), an inviscid CFD simulation module 
can predict the fluid properties for the inviscid outer region 
604 . 

FIGS. 6 a and 6b also depict a boundary -layer region 606 , 
located near a wing surface 602 . The boundary-layer region 
606 is typically characterized by a sharply increasing velocity 
profile 608 . Skin friction causes the velocity of the fluid flow 
very close to the wing surface to be essentially zero, with 
respect to the surface. A sharply increasing velocity profile 
608 develops as the velocity increases from a near-zero veloc- 
ity to the free-stream velocity. The sharply increasing velocity 
profile 608 in the boundary-layer region 606 creates shear 
forces within the botmdary -layer region 606 . Due to the inter- 
nal shear forces, viscous properties of the fluid influence fluid 
flow in the boundary-layer region 606 . Therefore, a simula- 
tion of the fluid flow in the boundary-layer region 606 must 
account for viscous contributions to the flow dynamics. In 
some cases, the fluid flow in boundary-layer region 606 may 
be characterized as turbulent flow (region 610 ). Due to fluid 
vorticity, viscous properties of the fluid influence the fluid 
flow. Thus, a simulation of the turbulent flow should also 
account for viscous contributions to the flow dynamics. For 
purposes of this discussion, laminar and turbulent regions are 
treated as one boundary-layer region and simulated using a 
single CFD simulation module. 

A CFD simulation module that accounts for viscosity may 
also be called a viscous CFD simulation module or a bound- 
ary-layer CFD simulation module. Exemplary equations for 
three-dimensional boundary-layer flow are: 
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kinetic energy thickness 0* is defined as: 


e* = f”(i 

Jo ' ' U c' 'PeUe 


where c p is the specific heat at constant pressure, p is the 
dynamic viscosity, k is the thermal conductivity, and T is the 
temperature. Often, the boundary-layer equations are further 111 
simplified, for instance by the sweep/taper theory where cer- 
tain derivatives in the spanwise direction can be neglected or 
approximated, allowing a three-dimensional boundary-layer 
solution to be obtained on a two-dimensional curved compu- 
tational mesh. Another simplification is to integrate out the 15 
details of the flow perpendicular to the surface of the airplane, 
resulting in an integral boundary-layer method, further 
described below with respect to equations 18 and 19. 

Below, exemplary field equations for a two-dimensional , f) 
integral boundary -layer CFD simulation module are provided 
according to a Drela boundary-layer technique. Drela, M. 
“XFOIL: An Analysis and Design System for Low Reynolds 
Number Airfoils,” pp. 1-12, Proceedings of the Conference 
on Low Reynolds Number Aerodynamics (T. J. Mueller ed., 25 
Univ. of Notre Dame, Notre Dame, Ind., 1989). 

Equation 18, below, represents a boundary- layer integral 
momentum equation for compressible flow: 
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Equation 18 


where 0 is the momentum thickness, H is the shape factor, M e , 
is the boundary-layer edge Mach number, u e is the boundary- 
layer edge velocity, and Cyis the skin friction coefficient. 

Equation 19, below, represents a boundary- layer kinetic 
energy integral equation: 

40 


dH * __ , e du e 

e— + (2H + H*(l - H))- — - = 2 C D 

CS X Vtg u X 


c, 

■»*f- 


Equation 19 


45 

As used in equations 18 and 19, above, shape factors H, H*, 
and H** are defined as: 



displacement thickness 8* is defined as: 



momentum thickness 0 is defined as: 


density thickness 8** is defined as: 


(5** = 


f< 


1 _ -dy; 
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skin friction coefficient C^is defined as: 



and dissipation coefficient C D is defined as: 

1 du 

C D = — 3 r— dy. 

P'<4 Jo dy 

Given the local velocity u, solving equations 1 8 and 1 9 for 
local displacement thickness 8* and momentum thickness 0, 
the boundary-layer CFD simulation module can predict the 
fluid properties for portions of the fluid flow within the 
boundary-layer region 606. This is the direct-mode simula- 
tion. Alternatively, the inverse-mode solution takes displace- 
ment thickness 8* as input and solves for local momentum 
thickness 0 and local velocity u. Additional characteristics of 
the boundary layer, including boundary -layer thickness, can 
also be determined once the fluid properties are known. 

In some cases, the boundary -layer CFD simulation module 
returns a transpiration flux value or transpiration velocity 
value representative of a boundary -layer region having a pre- 
dicted thickness. A transpiration flux is a fictitious fluid flow 
into or out of the computer-generated aircraft surface. The 
transpiration flux displaces a portion of the outer region 604 
in the same way as would the presence of a boundary layer 
having a predicted thickness. Using transpiration to represent 
the boundary layer may be more convenient. For instance, the 
inviscid CFD simulation module need only be modified 
where it calculates the surface boundary flow. It also allows 
the boundary-layer region to be discretized using a set of 
boundary-layer prediction points instead of using small vol- 
umes. In general, as the magnitude of the transpiration flux 
increases, the fictitious fluid flow increases, simulating a 
thickening boundary layer. In some cases, the transpiration 
flux can be used to create a fictitious flow of air into the 
aircraft surface (negative flux), thereby simulating a bound- 
ary layer having a reducing thickness. 

The transpiration flux can be determined using the output 
of the Drela boundary -layer teclmique described in equations 
1 8 and 19, above. For example, the transpiration flow velocity 
v w of the transpiration flux can be determined using: 


9 = 
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Equation 20 
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where p e is the density of the fluid flow at the boundary-layer 
edge, u e is the velocity of the fluid flow at boundary-layer 
edge, 8* is the computed boundary-layer displacement thick- 
ness, and s is an arc length along the computer-generated 
aircraft surface. In some cases, the arc length s is the average 
distance between boundary-layer prediction points on the 
computer-generated aircraft surface. Equation 20 is taken 
from Lock, R. C., and Williams, B. R., “Viscous-Inviscid 
Interactions in External Aerodynamics,” Prog. Aerospace 
Sci., Vol. 24, 1987, pp. 51-171. Thus, the transpiration mass 
flux (density p e times the transpiration flow velocity v w ) is 
equal to the rate of change of the product of the local density 
p e , local velocity u e , and boundary-layer displacement thick- 
ness 8* along the solution strip. A finite difference method 
can be used to compute the derivative in equation 20. For 
example, neighboring boundary-layer prediction points 
along the surface of the computer-generated aircraft surface 
can be used with a second order, backward Lagrange poly- 
nomial formulation to compute the derivative values. 

For the case of three-dimensional or sweep/taper bound- 
ary-layer flow solutions, there is an analog to equation 20 that 
has two terms on the right-hand side. The first term with the 
streamwise displacement thickness and the second with the 
crossflow displacement thickness. 

The results of the inviscid CFD simulation module (using, 
for example, equation 13) are used by a viscous or boundary- 
layer CFD module (using, for example equations 18, 19, and 
20) to predict the boundary-layer fluid properties for each 
boundary-layer prediction point. The boundary-layer fluid 
properties (e.g., boundary-layer thickness or transpiration 
flux) can then be used by the inviscid CFD simulation module 
to predict an updated or refined set of inviscid fluid property 
values. In this way, the inviscid and viscous CFD simulation 
modules can be used to iterate or refine the results of the 
computer-generated simulation. 

5. Computer and Computer Network System 

The embodiments described herein are typically imple- 
mented as computer software (computer-executable instruc- 
tions) executed on a processor of a computer system. FIG. 7 
depicts an exemplary computer system 700 configured to 
perform any one of the above-described processes. Computer 
system 700 may include the following hardware components: 
processor 702, data input devices (e.g., keyboard, mouse, 
keypad) 704, data output devices (e.g., network connection, 
data cable) 706, and user display (e.g., display monitor) 708. 
The computer system also includes non-transitory memory 
components including random access memory (RAM) 710, 
hard drive storage 712, and other computer-readable storage 
media 714. 

Processor 702 is a computer processor capable of receiving 
and executing computer-executable instructions for perform- 
ing any of the processes described above. Computer system 
700 may include more than one processor for performing the 
processes. The computer-executable instructions may be 
stored on one or more types of non-transitory storage media 
including RAM 710, hard drive storage 712, or other com- 
puter-readable storage media 714. Other computer-readable 
storage media 714 include, for example, CD-ROM, DVD, 
magnetic tape storage, magnetic disk storage, solid-state stor- 
age, and the like. 

FIG. 8 depicts an exemplary computer network for distrib- 
uting the processes described above to multiple computers at 
remote locations. One or more servers 810 may be used to 
perform portions of the process described above. For 
example, one or more servers 810 may store and execute 
computer-executable instructions for receiving information 
for generating a computer-generated simulation. Hie one or 
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more servers 810 are specially adapted computer systems that 
are able receive input from multiple users in accordance with 
a web-based interface. Hie one or more servers 810 are able 
to communicate directly with one another using a computer 
5 network 820 including a local area network (LAN) or a wide 
area network (WAN), such as the Internet. 

One or more client computer systems 840 provide an inter- 
face to one or more system users. The client computer sys- 
tems 840 are capable of communicating with the one or more 
to servers 810 over the computer network 820. In some embodi- 
ments, the client computer systems 840 are capable of run- 
ning a web browser that interfaces with a web-enabled system 
running on one or mover server machines 810. The web 
browser is used for accepting input data from the user and 
15 presenting a display to the user in accordance with the exem- 
plary user interface described above. The client computer 840 
includes a computer monitor or other display device for pre- 
senting information to the user. Typically, the client computer 
840 is a computer system in accordance with the computer 
20 system 700 depicted in FIG. 7. 

Although the invention has been described in considerable 
detail with reference to certain embodiments thereof, other 
embodiments are possible, as will be understood by those 
skilled in the art. 

25 We claim: 

1 . A computer-implemented method of generating a fluid- 
flow simulation over a computer-generated surface using one 
or more processors, the simulation including an inviscid 
fluid-flow mesh comprised of a plurality of inviscid fluid cells 

30 and a viscous fluid-flow mesh comprised of a plurality of 
boundary-layer fluid cells, at least some of the boundary- 
layer fluid cells being on or adjacent to the computer-gener- 
ated surface, the method comprising: 

determining, using the one or more processors, an initial 
35 fluid property, for at least one inviscid fluid cell using an 

inviscid fluid simulation that does not simulate fluid 
viscous effects; 

determining, using the one or more processors, an initial 
boundary-layer fluid property for at least one of the 
40 boundary-layer fluid cells using the initial fluid property 
and a viscous fluid simulation that simulates fluid vis- 
cous effects, 

wherein the at least one inviscid fluid cell is located in 
relation to the at least one boundary-layer fluid cell 
45 such that anupdatedboundary-layerfluidproperty for 

the at least one boundary-layer fluid cell is influenced 
by the at least one inviscid fluid cell; 

determining, using the one or more processors, the updated 
boundary-layer fluid property for the at least one bound- 
50 ary-layer fluid cell using the initial fluid property, initial 
boundary-layer fluid property, and an interaction law, 
wherein the interaction law approximates the inviscid 
fluid simulation using a matrix of aerodynamic influ- 
ence coefficients computed using a two-dimensional 
55 surface panel technique and a fluid-property vector; 

determining, using the one or more processors, an updated 
fluid property for the at least one inviscid fluid cell using 
the updated boundary-layer fluid property. 

2. The computer-implemented method of claim 1, wherein 
60 the initial fluid property is an initial fluid velocity and the 

updated fluid property is an updated fluid velocity. 

3. The computer-implemented method of claim 1, wherein 
the initial boundary-layer fluid property is an initial displace- 
ment thickness and the updated boundary-layer fluid property 

65 is an updated displacement thickness. 

4. The computer-implemented method of claim 1, wherein 
the matrix of aerodynamic influence coefficients is adapted to 
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account for compressibility by dividing each row by (1- 
Mj 2 ) 172 , wherein Mj is the local Mach number. 

5. The computer-implemented method of claim 1, wherein 
a fluid property of the fluid property vector is determined 
using the following relation to simulate supersonic flow con- 
ditions: 


u e(n+l) = ” 


1 l dm M 
ds 


<Jm (n+ 1, -j 

ds y 


where u e(M+1) is the edge velocity at time step n+1, u e(n) is 
the edge velocity at time step n, M e , is the local Mach 
number, and m, , and m (n+1) is the product of: edge 
density p e , edge velocity u e , and boundary-layer thick- 
ness 6* at time step n and time step n+1, respectively. 
6. The computer-implemented method of claim 1, wherein 
a fluid property of the fluid property vector is determined 
using the following relation to simulate supersonic flow con- 
ditions over a swept/tapered wing: 

1 

w e(/i+l) — M e(n) . =( V *<«+1) — 

V Ml - 1 


where u e(H+/) is the edge velocity at time step n+1, u e(M) is 
the edge velocity at time step n, M ei is the local Mach 
number, and v M , (H+1) and v u , (w) are the transpiration 
velocities at the surface at time step n and time step n+1 , 
respectively. 

7. The computer-implemented method of claim 1, wherein 
the matrix of aerodynamic influence coefficients is adapted 
row-by-row to implement either the subsonic or supersonic 
interaction law depending on the local Mach number. 

8. The computer-implemented method of claim 7, wherein 
the matrix of aerodynamic influence coefficients is further 
adapted to utilize a linear weighted average of the subsonic 
and supersonic interaction laws to approximate the changes 
to the inviscid fluid flow for local Mach numbers near Mach 
1. 

9. The computer- implemented method of claim 8, wherein 
the linear weighted average is used for local Mach numbers 
ranging between about 0.97 and about 1 .03. 

10. The computer-implemented method of claim 1, 
wherein the matrix of aerodynamic influence coefficients is 
adapted using a curvature correction to approximate a bound- 
ary layer characterized as having a thickness and a curvature. 

11. The computer-implemented method of claim 1, 
wherein the interaction law approximates the inviscid fluid 
simulation using a matrix of aerodynamic influence coeffi- 
cients computed using an axisymmetric surface panel tech- 
nique instead of a two-dimensional surface panel technique. 

12. The computer-implemented method of claim 1, 
wherein the computer-generated surface is a computer-gen- 
erated wing surface of an aircraft. 

13. A non-transitory computer-readable storage medium 
comprising computer-executable instructions for generating 
a fluid-flow simulation over a computer-generated surface, 
the simulation including an inviscid fluid-flow mesh com- 
prised of a plurality of inviscid fluid cells and a viscous 
fluid-flow mesh comprised of a plurality of boundary-layer 
fluid cells, at least some of the boundary-layer fluid cells 
being on or adjacent to the computer-generated surface, the 
instructions for: 


determining an initial fluid property, for at least one invis- 
cid fluid cell using an inviscid fluid simulation that does 
not simulate fluid viscous effects; 
determining an initial boundary-layer fluid property for at 
5 least one of the boundary-layer fluid cells using the 
initial fluid property and a viscous fluid simulation that 
simulates fluid viscous effects, 

wherein the at least one inviscid fluid cell is located in 
relation to the at least one boundary-layer fluid cell 
such that an updated boundary-layer fluid property for 
the at least one boundary-layer fluid cell is influenced 
by the at least one inviscid fluid cell; 
determining the updated boundary-layer fluid property for 
j 5 the at least one boundary-layer fluid cell using the initial 

fluid property, initial boundary-layer fluid property, and 
an interaction law, 

wherein the interaction law approximates the inviscid 
fluid simulation using a matrix of aerodynamic influ- 
20 ence coefficients computed using a two-dimensional 

surface panel technique and a fluid-property vector; 
determining an updated fluid property for the at least one 
inviscid fluid cell using the updated boundary-layer fluid 
property. 

25 14. The computer-readable medium of claim 13, wherein 

the initial fluid property is an initial fluid velocity and the 
updated fluid property is an updated fluid velocity. 

15. The computer-readable medium of claim 13, wherein 
the initial boundary -layer fluid property is an initial displace- 

30 ment thickness and the updated boundary -layer fluid property 
is an updated displacement thickness. 

16. The computer-readable medium of claim 13, wherein 
the matrix of aerodynamic influence coefficients is adapted to 
account for compressibility by dividing each row by (1- 

35 Mj 2 ) 1 2 , wherein M ; is the local Mach number. 

17. The computer-readable medium of claim 13, wherein a 
fluid property of the fluid property vector is determined using 
the following relation to simulate supersonic flow conditions: 


40 
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45 where u e( „ +1) is the edge velocity at time step n+ 1 , u e(M) is the 
edge velocity at time step n, M ei is the local Mach number, and 
m ( , ;) and m( n+1) is the product of: edge density p e , edge veloc- 
ity u e , and boundary-layer thickness 8 * at time step n and time 
step n+1 , respectively. 

50 18. The computer-readable medium of claim 13, wherein a 

fluid property of the fluid property vector is determined using 
the following relation to simulate supersonic flow conditions 
over a swept/tapered wing: 


u e(n+ 1) - U <(n, , “ l Vr;T 

Ml- 1 


60 where u e( „ +1 j is the edge velocity at time step n+1, u e(n) is the 
edge velocity at time step n, M e! . is the local Mach number, and 
v >, <n+i > and v w( „, are the transpiration velocities at the surface 
at time step n and time step n+1 , respectively. 

19. The computer-readable medium of claim 13, wherein 
65 the matrix of aerodynamic influence coefficients is adapted 
row-by-row to implement either the subsonic or supersonic 
interaction law depending on the local Mach number. 
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20. The computer-readable medium of claim 19, wherein 
the matrix of aerodynamic influence coefficients is further 
adapted to utilize a linear weighted average of the subsonic 
and supersonic interaction laws to approximate the changes 

to the inviscid fluid flow for local Mach numbers near Mach 5 

1 . 

21. The computer-readable medium of claim 20, wherein 
the linear weighted average is used for local Mach numbers 
ranging between about 0.97 and about 1 .03. 

22. The computer-readable medium of claim 13, wherein to 
the matrix of aerodynamic influence coefficients is adapted 
using a curvature correction to approximate a boundary layer 
characterized as having a thickness and a curvature. 

23. The computer-readable medium of claim 13, wherein 
the interaction law approximates the inviscid fluid simulation 15 
using a matrix of aerodynamic influence coefficients com- 
puted using an axisymmetric surface panel technique instead 
of a two-dimensional surface panel technique. 

24. The computer-readable medium of claim 13, wherein 
the computer-generated surface is a computer-generated 20 
wing surface of an aircraft. 



